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We address the energetic stability of the graphene/SiC(0001) interface and the associated binding 
mechanism by studying a series of low-strain commensurate interface structures within a density 
functional scheme. Among the structures with negligible strain, the 6\/3x6\/3 -R30 SiC periodicity 
shows the lowest interface energy, providing a rationale for its frequent experimental observation. 
The interface stability is driven by the enhanced local reactivity of the substrate-bonded graphene 
atoms undergoing sp 2 -to-sp 3 rehybridization (pyramidalization) . By this mechanism, relaxed struc- 
tures of higher stability exhibit more pronounced graphene corrugations at the atomic scale. 
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Epitaxially grown graphene on SiC provides a promis- 
ing solution for the realization of high-performance 
carbon-based electronic devicesP^On the Si-face of SiC, 
the interface between the epitaxial graphene overlayers 
and the SiC substrate consists of a thin carbon buffer 
l a y ei EHll which has bee n char acterized through several 
experimental techniques! 1 * 3 * 6 * 9 * *^! Angle-resolved photoe- 
mission spectroscopy data show the presence of well- 
developed a bands and are thus consistent with a buffer 
layer having the same topology as graphene! 9 ^ However, 
the absence of the free-standing graphene Dirac cones in 
the 7T bands points to the occurrence of covalent bond- 
ing between the buffer layer and SiC (000 1) j^ l as a l so 
corroborated by x-ray photoemission dataPESl A moire 
pattern of 6V 3 x 6-^/3 i?30° SiC periodicity usually arises 
in electron d iffrac tio n 1 * 9 " 11 * and in scanning tunneling 
measurements^^ 3 -* due to the lattice constant mismatch 
between graphene and SiC. The observed moire is com- 
monly understood as a coincidence lattice between a 
6V3x6\/3 -R30 Si^surface reconstruction and a 13 x 13 
graphene sheet.* 6 ^* 14 * 15 * However, several other SiC pe- 
riodicities, such as 2\fZ x 2v / 3i?30°, 4x4, 5x5, 
are compatible with an almost strain-free graphene/SiC 
interfacef^iHIIl thoug h only few of them have been found 
in experiments! 11 * 18 * A further indication that the exis- 
tence of a low-strain coincidence lattice is not a suffi- 
cient condition for its experimental realization can be 
inferred from the observation of a 2.2°-rotated commen- 
surate graphene phase on the C facef^l which does not 
have a counterpart on the Si face. Hence, legitimate ques- 
tions arise about the underlying physical mechanisms 
which drive the occurrence of specific structures at the 
graphene/SiC(0001) interface. 

In this Rapid Communication, we address the energet- 
ics of the graphene/SiC(0001) interface by considering a 
series of model structures differing by the SiC periodic- 
ity, the strain in the graphene layer, and the rotation 
angle between the graphene and SiC lattices. Our study 
shows that the experimentally observed 6-\/3 x 6-\/3 i?30° 
periodicity is favored over other low-strain commensu- 
rate structures by its higher binding energy. An analysis 



of the binding energies in terms of C-Si covalent bond 
energies reveals that the driving mechanism leading to 
interface stability is dominated by the sp 2 -to-sp 3 rehy- 
bridization undergone by the substrate-bonded graphene 
atoms. As a consequence, a clear connection is estab- 
lished between the binding energy and the atomic-scale 
corrugations of the graphene layer. 

In our density functional theory calculations, the 
exchange-correlation energy is described within a gener- 
alized gradient approximation (GGA)|251 We also checked 
that the description of va n der Waals interactions be- 
yond the GGA leve l 21 * 22 -* does not change our conclu- 
sions. The valence wave functions and the electron charge 
are expanded in plane wave basis sets defined by ki- 
netic energy cutoffs of 35 and 350 Ryp^ respectively. 
The core-valence interaction is described through ultra- 
soft pseudopotentialsP^l Our slab geometries comprise 
one graphene layer and five bilayers of SiC, of which the 
bottom C atoms are saturated with H atoms. The slab 
is separated from its periodic replicas by vacuum regions 
of 30 A. For each interface system the sampling of the 
Brillouin zone is at least equivalent to a 10 x 10 fc-point 
mesh in the primitive cell of graphene, which ensures er- 
rors lower than 0.01 eV on the total energy per surface 
Si atom. All atoms but those in the bottom two SiC 
bilayers are allowed to relax. We used the Quantum 
ESPRESSO packageP 

Focusing on to x m and mVZxmVsmO SiC period- 
icities (denoted as Nto and Rto, respectively), we search 
for commensurate interface structures of tractable size in 
which the graphene sheet is subject to limited strain, ei- 
ther tensile or compressive.^ This procedure results in 
a series of interface structures, each of them character- 
ized by an integer m and a rotation angle between the 
graphene and SiC lattices (Table [I]) . Structures with the 
same SiC periodicity but with different rotation angle 
are distinguished by primes. The largest SiC period- 
icity considered here corresponds to the experimentally 
observed 6-\/3 x 6\/3 i?30°, requiring up to 1526 atoms 
in the simulation cell. For all retained structures, the 
graphene strain is below 1.5% (cf. Table |l| except for Rl 
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(8.3%). The latter presents an unrealistic strain level, 
but has been included here for comparison with previous 
work!^ For each structure, we also report in Table [i] the 
calculated interface energy per surface Si atom E- lu x, re- 
ferred to ideal, unstrained graphene and unreconstructed 
SiC(OOOl). Apart from the Rl interface, all interface 
structures are energetically stable (E- lnt < 0). 

As shown in Fig. [TJ the energies of interfaces with 
strain lower than 0.2% decrease monotonically with in- 
creasing graphene/SiC rotation angle. The lowest energy 
is found for the 6-\/3x6-\/3 i?30° periodicity (R6), provid- 
ing an energetic argument for its frequent experimental 
observation. Semi-empirical corrections accounting for 
the van der Waals interactions^ just shift the interface 
energies by a constant amount and thus do not affect this 
comparison. With respect to the low-strain interfaces 
with close rotation angles, the structures with higher 
tensile and compressive strain yield lower and higher in- 
terface energies, respectively. This behavior is consis- 
tent with the dependence of the chemical reactivity of 
graphene on strain,^ leading in particular to a very sta- 
ble R4" interface. However, to the best of our knowledge, 
this periodicity has not been observed experimentally, 
suggesting that such levels of graphene strain cannot be 
sustained during growth. 

To investigate the factors playing a role in the interface 
stability, the interface energy has been decomposed as 
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where E str is the energy cost to strain ideal graphene 
without deformations, i?def is the energy needed to de- 
form the strained graphene layer into the final atomic 
configuration, and -E s tick is the energy gained by sticking 



TABLE I: Energetics and structural properties of the re- 
laxed graphene/SiC(0001) interfaces: graphene/SiC crystal- 
lographic rotation angle a, strain s on the graphene sheet, 
interface energy per surface Si atom Ei nt , and vertical spread 
of the graphene C atoms Azg- E- m t is decomposed in strain 
energy E s ti, deformation energy Edei, and sticking energy 
Sstick of the graphene layer. The interfaces are separated 
in three groups according to strain: \s\ ^ 0.2%, s > 0.2%, 
and s < —0.2%. Energies are in eV and Azg in A. 
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FIG. 1: (Color online) Interface energy E[ nt vs graphene/SiC 
rotation angle a. The structures are distinguished by the 
strain s in the graphene layer: very weak strain with \s\ ^ 
0.2% (circles), tensile strain with s > 0.2% (diamonds), and 
compressive strain with s < —0.2% (squares). The dashed 
line is a linear fit of Ei nt for the structures with Isl ^ 0.2%. 



the so-obtained graphene layer onto SiC. In this decom- 
position, the deformation energy of the SiC substrate is 
included in -Estick and is not singled out because its con- 
tribution is found to be negligible (< 0.02 eV per surface 
Si atom) . The calculated values of the three components 
of the binding energy in Eq. ([!]) are given in Table [I] 
Focusing first on interfaces with strain lower than 0.2%, 
we notice that the strain energy is much smaller than the 
interface energy and thus does not influence the stability. 
Furthermore, one observes that Ei nt results from a com- 
pensation effect between the cost of -Edef and the gain 
of -Estick, with S s tick roughly scaling with E dc{ . In turn, 
this directly connects the interface energy with the de- 
formation energy, resulting in higher stabilities for larger- 
degrees of graphene deformation. Overall, these observa- 
tions qualitatively also hold for the interfaces with larger 
strains except for Rl, in which the considerable strain 
energy makes the structure unstable. 

To better understand their relation with the inter- 
face energy, we analyze in more detail the nature of 
the deformations undergone by graphene. It has been 
suggested that the structural properties of the interface 
could be driven by dangling bond saturation of the sur- 
face Si a toms through bond formation with the graphene 
C atomsPED Since the ideal moire pattern also presents 
some surface Si ato ms which do not have any C atom 
directly atopJM-iEGD one CO uld then argue that deforma- 
tions are mainly induced by lateral adjustments of the 
C atoms. However, it has been shown that the interface 
energy as described in a bond energy picture does not un- 
dergo significant variations when comparing the case of 
undistorted graphene with that of the fully relaxed buffer 
layerP^ Furthermore, when focusing on interface struc- 
tures with \s\ < 0.2% (Table [T]), one observes that the 
deformation energy approximately scales with the verti- 
cal spread of the graphene C atoms AzgW^ suggesting 
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FIG. 2: (Color online) Distribution of the pyramidalization 
angle P of the C atoms in the buffer layer for the inter- 
face structures with s ^ 0.2%. The solid bars refer to Si- 
bonded graphene C atoms identified by a cutoff distance of 
2.3 A, while the open bars refer to non-bonded C atoms. The 
graphene/SiC rotation angle a, the interface energy Ei nt , and 
the fraction /si of Si-bonded graphene C atoms are indicated 
for each structure. In the bottom panel, the inset illustrates 
the 7r-axial vector v-n- and the pyramidalization angle 9 P . 



atom with sp 2 hybridization as in the ideal honeycomb 
structure, and takes the value of 19.47° for a C atom 
in the tetrahedral sp 3 configuration. The distributions 
of the calculated pyramidalization angles of the C atoms 
in the buffer layer are shown in Fig. [2] for the interface 
structures with \s\ < 0.2%. For all the structures consid- 
ered here, the angle distribution is bimodal. To elucidate 
this property, we separate the graphene C atoms in two 
groups by identifying those forming bonds to surface Si 
atoms (Fig. [2]). This decomposition shows that graphene 
C atoms forming a C-Si bond are responsible for the peak 
at high pyramidalization angles, while the non-bonded C 
atoms give rise to the feature at lower 9 p . 

The fraction of Si-bonded graphene atoms does not 
vary significantly among the considered interface struc- 
tures (Fig. [2| and cannot explain the differences in the 
calculated interface energies, as noted above. In this per- 
spective, it is important to observe that the peak at high 
pyramidalization angles shifts significantly toward higher 
values as interfaces of increasing stability are considered. 
This suggests a direct link between the local C-Si bond 
energy and the pyramidalization angle of the involved 
C atom. Indeed, the bond energy of an individual C-Si 
bond can be expr essed in terms of the rehybridization of 
the C 7r-orbital asPl2l 
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where e gr is the energy cost for isolating a single it orbital 
from the conjugated tt system of graphene, and the ma- 
trix elements e s and e p involve surface Si orbitals and C 
orbitals of s andp z character, respectively. Since e s < e p , 
C atoms with larger 9 P are more reactive and lead to 
stronger C-Si bonds. 

To evaluate the bond energy variations as a function 
of dp , we use an atomistic model consisting of a pyrami- 
dal structure in which a single surface Si atom is exposed 
to the graphene layer as illustrated in Fig. |3ja). In 
the range of interest (9 p < 20°), the dependence on 9 P is 
found to be approximately linear [Fig. [3jb)] , as expected 



that the main contributions come from out-of-plane de- 
formations rather than from lateral displacements. 

In view of these considerations, we further analyze the 
deformations by focusing on the rehybridization under- 
gone by the graphene C atoms. Rehybridization effects 
have already been pointed out to explain the pinning 
of graphene to a metallic substrate upon deposition of 
adsorbatesp3 To this purpose, we consider the pyrami- 
dalization angle 9 p (Ref. [55]) formed by each of the C 
atoms in the buffer layer. As shown in the inset in Fig. 
[2j the definition of 9 P is based on the construction of the 
7r-orbital axis vector vv, which forms equal angles 9 a7r 
with the three cr-bonds connecting the central C atom 
to its graphene neighbors. The pyramidalization angle 
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FIG. 3: (Color online) (a) Atomistic configuration used to 
compute -Ec-Si as a function of the pyramidalization angle 
9 P [only the top of the SiC pyramid construction (Ref. I16[) 
is shown], (b) Calculated -Ec-Si vs 9 P and fit according to 
Eq. S (dashed line). We used a fixed C-Si distance de-Si of 
2.05 A, corresponding to typical C-Si bond distances at the 
graphene/SiC(0001) interface. 
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FIG. 4: (Color online) Sticking energy -E st i c k of the graphene 
buffer layer vs model energy -E ro h y b accounting for the C tt- 
orbital local rehybridizations. Both energies are given per 
surface Si atom. The structures are distinguished according 
to strain, as in Fig. [TJ 

from Eq. The calculated Bc-Si values are well re- 
produced by Eq. ([2| with e gr = 2.8 eV, e s = — 5.6eV and 
e p = — 1.9 eV. In this way, we obtain for each Si-bonded 
graphene C atom a bond energy Ec_si(#p) from its 9 p 
value in the relaxed interface structure. This allows us 
to define a model energy -Erehyb based on the sum of all 
these individual contributions which account for the local 
trigonal hybridization of the C ir orbitals. 

In Fig. |4| we confront -Erehyb with the sticking energy 
E s tick for all structures in Table [j] but RlP3 The struc- 



tures with graphene strain \s\ < 0.2% show that the vari- 
ations in E s tick are fully captured by the variations of 
Erehyb- This correlation approximately also holds when 
considering the structures of higher strain (|s| > 0.2%), 
indicating that the C reactivity is affected by strain ef- 
fects to a lesser extent. These considerations demonstrate 
that the interfacial binding mechanism is directly related 
to the atomic-scale corrugation of the graphene and to 
the associated enhancement of the chemical reactivity of 
the Si-bonded graphene C atoms due to their partial sp 2 - 
to-sp 3 rehybridization. 

In conclusion, we highlight the role of carbon rehy- 
bridization in the binding of the graphene layer to the 
SiC(0001) substrate. This notion carries general valid- 
ity and is expected to hold for other graphene interfaces 
involving covalent bonding. The pyramidalization of the 
carbon atoms also connects with the atomic-scale corru- 
gation of the graphene layer and could further lead to 
a spatial modulation of the graphene reactivity toward 
adsorbates. 
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